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Geometrically frustrated magnetic molecules have attracted a lot of interest in the field of molecular mag- 
netism as well as frustrated Heisenberg antiferromagnets. In this article we demonstrate how an approximate 
diagonalization scheme can be used in order to obtain thermodynamic and spectroscopic information about 
frustrated magnetic molecules. To this end we theoretically investigate an antiferromagnetically coupled spin 
system with cuboctahedral structure modeled by an isotropic Heisenberg Hamiltonian. 
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1. Introduction 

The complete understanding of small magnetic systems such as magnetic molecules is compulso- 
rily connected to the knowledge of their energy spectra. From the energy spectra all spectroscopic, 
dynamic, and thermodynamic properties of the spin systems can be obtained. Unfortunately, an 
exact calculation of the spectrum is often restricted due to the huge dimension of the Hilbcrt space 
even if one works within the most simple isotropic Heisenberg model. The dimension grows for a 
system of N spins with spin quantum number s exponentially and is (2s + \) N . 

In order to get insight into the properties of large magnetic molecules one can access several 
numerical methods which have been developed in the past. Of course, the ultimate method of choice 
would be an exact numerical diagonalization yielding the complete energy spectrum. In recent years 
there has been enormous progress on extending the range of applicability of the exact numerical 
diagonalization of the Heisenberg model. To this end the use of spin-rotational symmetry [[TJ[2] in 
combination with point-group symmetries [ [3l [4] [6] can be of great advantage with respect to a 
reduction of computational requirements, i.e. a need of hardware resources and computation time. 
Apart from the exact numerical diagonalization technique the magnetism of magnetic molecules can 
be very well investigated using complementary methods such as Density Matrix Renormalization 
Group (DMRG) [E1IH1I5], Lanczos [[TO], or Quantum Monte Carlo (QMC) [QI1HQ31 techniques. 
Nevertheless, also these methods suffer from theoretical limitations, QMC for instance in systems 
with geometric frustration. 

Currently magnetic molecules which exhibit geometric frustration are of special interest due to 
the richness of physical phenomena like plateaus and jumps of the magnetization for varying field as 
well as special features of their spectra such as low- lying singlets [[T21[T51[TEl[r7] . In this respect a lot 
of insight has been obtained by investigating molecular representations of archimedean-type spin 
systems [ 1181 119) . i.e. systems in which participating spins occupy the vertices of an archimcdean 
solid. Such representations have already been synthesized several years ago and exist for example as 
{Cu 12 La 8 } [[101 (cuboctahedron, s = 1/2) and {Mo 72 Fe 3 n} [EH1I221 (icosidodecahedron, s = 5/2). 
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In this context the molecular compound {Mo72Fc3o} is probably one of the most investigated 
magnetic molecules. However, a theoretical explanation of its special physical properties was so far 
given mostly by considering purely classical models [ 1231 I24j or directly related quantum mechanical 
counterparts like the rotational-band model [ [25] [26] . 

In this paper we want to show how the approximate diagonalization technique which has been 
developed and applied to unfrustratcd, i.e. bipartite, magnetic molecules in Ref. [ |S] can be used 
in order to determine the energy spectrum of geometrically frustrated magnetic molecules. The 
idea of this technique is to diagonalize the full Hamiltonian in a reduced basis set. The basis set 
itself is an eigenbasis of the rotational-band Hamiltonian. Such an ansatz was already used by 
Oliver Waldmann [[27] in order to interpret inelastic neutron spectra of {Mo72Fe3o} [ 118]. We 
will demonstrate that in contrast to bipartite systems for a frustrated spin systems the rotational- 
band states of all non-trivially different sublattice colorings have to be taken into account in 
order to achieve a reliable convergence of energy levels. Throughout this paper we use a spin 
system with cuboctahedral structure and spin quantum numbers s = 1 and s = 3/2 as well as an 
antiferromagnetic coupling as an archetypical example of a frustrated magnetic spin system. 

This paper is organized as follows. In Sec. [5] the general theoretical description of the system 
within the Hciscnberg model as well as general remarks on the use of irreducible tensor opera- 
tors and point-group symmetries are given. In Sec. [3] the theoretical basics of the approximate 
diagonalization within the isotropic Hciscnberg model are briefly reviewed and specified for the 
cuboctahedral spin system. The convergence behaviour of the approximate diagonalization is dis- 
played and deeply discussed as well as the specific heat and zero-field magnetic susceptibility for a 
cuboctahedron with s = 3/2. Furthermore, a numerically based finding of an approximate selection 
rule is reported. This paper closes with a Summary in Sec. [4] 



2. Theoretical method 

In order to model the physics of antiferromagnetic molecules it has been shown that an isotropic 
Hciscnberg Hamiltonian with an additional Zccman-term and an antiferromagnetic nearest-neighbor 
coupling provides the dominant terms. Such a Hamiltonian looks like 

# = - E J «£W ■ SC?) + 5MB§ -B. (1) 

<i,i> 

The indices of the sum are running over all pairs < i, j > of interacting spins i and j. The first 
part consisting of the sum over single spin operators s(i) at sites i interacting with the coupling 
strength Jy < refers to the Heisenberg exchange whereas the second part - the Zeeman-term - 
couples the total spin S to an external magnetic field B. 

Taking without loss of generality = J for interacting spins the Heisenberg part assumes the 
following form 

^Hciscnberg = ~J ^ ' £(•?')> ( 2 ) 

<i,i> 

where each coupling is counted only once. Since due to SU(2) symmetry the commutation relations 
^Hciscnberg, S =0 hold, a common eigenbasis of i?Heisenberg, S 2 and S z can be found and 

the effect of an external magnetic field B = B ■ e z can be included later, i.e. 

E v {B) = E„ + g(j, B BM v . (3) 

Here E v denotes the energy eigenvalues, \v) the eigenstates, and M v denotes the corresponding 
magnetic quantum number. 

For the matrix representation of the Heisenberg Hamiltonian ^ the irreducible tensor operator 
technique is used [[TJ[2][6]. The Heisenberg Hamiltonian is expressed in terms of irreducible tensor 
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operators and its matrix elements are evaluated using the Wigner-Eckart-theorem 
(aSM\T^\a' S' M') = 

(-l) s - M (aS\\TM\\a>S>)( _^ J ^ ) . (4) 

Equation ((4]) states that a matrix clement of the g-th component of an irreducible tensor operator 
T( fc ) of rank k is given by the reduced matrix element {a S\ \ \a' S') and a factor containing a 
Wigner-3J symbol [HE]- The basis in which the Hamilton matrix is set up is of the form \aS M). 
a refers to a set of intermediate quantum numbers given by addition rules when coupling single 
spins s(i) to a total spin S with spin quantum number S. The underlying spin-coupling scheme 

directly influences the form of the irreducible tensor operator and further on the successive 
calculation process of the reduced matrix elements in Eq. (j4|). 

By using irreducible tensor operators and the Wigner-Eckart-theorem it is possible to drastically 
reduce the dimensionality of the problem, i.e. of the Hamilton matrices which have to be diago- 
nalizcd numerically. The total Hilbert space TC can be decomposed into subspaccs H(S,M = S). 
Such a decomposition results in block factorizing the Hamilton matrix where each block can be 
labelled by the total spin quantum number S. 

Additionally point-group symmetries lead to a further reduction of the matrices. Apart from 
that, states are labelled by the irreducible representation they are belonging to. This is very 
helpful in understanding the physics of the system. The incorporation of these symmetries results 
in symmetrized basis states which are constructed by the projection operator [ [29] 

V^\aSM)= (|E(x (n) (^))* Q( R H \ aSM ) ■ ( 5 ) 

Here l n denotes the dimension of the ro-th irreducible representation of the point-group Q which 
is of order h. G(R) refers to the symmetry operations of Q and x (R) denotes its character with 
respect to n. The effect of the symmetry operation G{R) on basis states of the form \aSM) is 
discussed in Refs. [03 El [6]. 

3. Approximate diagonalization for frustrated systems: the cuboctahedron 



Figure 1. Geometrical structure of the cuboctahedron [130]. 

In this section we follow a method of calculating approximate cigenstates and eigenvalues of 
the Hciscnbcrg Hamiltonian in Eq. ([2|) and determine approximately the energy spectrum of a 
spin system with cuboctahedral structure. In this system 12 spins of spin quantum number s 
occupy the vertices of a cuboctahedron interacting along its edges. The geometrical structure 
of the cuboctahedron is shown in Fig. [TJ A detailed description of the approximation scheme is 
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given in Rcf. [ 6] where it was successfully applied to determine approximate energy spectra and 
thermodynamic properties of spin rings, i.e. bipartite systems. This approximation rests on the 
idea to diagonalizc the full Hamiltonian in a reduced basis set. The basis set itself is an eigenbasis 
of the rotational-band Hamiltonian which from the point of view of perturbation theory can be 
understood as an approximation to the full Hamiltonian. For bipartite systems this (zeroth order) 
approximation is already very good [ [25] . 

The Heiscnbcrg Hamiltonian can be decomposed into two parts like 

H = Hrb + H' , (6) 

where -Hrb is the rotational-band Hamiltonian [ OH] and H' is an operator containing 

the remaining terms. The rotational-band Hamiltonian which is an effective quantum mechanical 
Hamiltonian based on classical assumptions looks like 

(7) 

Here N denotes the number of spins within the system and N s the number of sublattices the 
classical ground state of the system is composed of. The prefactor —DJ/(2N) can be seen as the 
effective coupling strength between the total spin S and the sublattice spins S„ which arise from 
coupling all single spins s(i) belonging to the n-th sublattice. The value of D directly depends on 
the system. It is chosen such that the energy of the ferromagnetic state of the system is matched; 
for the cuboctahedron it is D = 6. 

The full Heiscnbcrg Hamiltonian is now diagonalizcd within a reduced set i = 1, • • ■ , n re d, 

of basis states of H^b yielding approximate eigenstates and eigenvalues of H . The set of approxi- 
mate basis states is energetically ordered. 
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Figure 2. Part of the energy spectra of the rotational-band Hamiltonian for a cuboctahedron 
with twelve spins s = 1 (left) and s — 3/2 (right). Seven super-bands are colored. 



Before discussing the results of our approximate diagonalization we would like to characterize 
the eigenbasis of -Hrb • Figure [2] shows the low-lying part of the energy spectra of -Hrb for a 
cuboctahedron with s = 1 and s = 3/2 given by the rotational-band model. They exclusively 
consist of parabolas - so-called rotational bands. The eigenvalues Erb (Si , S2, S3 , S) of H rb depend 
only on the spin quantum numbers of the sublattice spins S n , with n = 1,2,3, and of the total spin 
S. Corresponding eigenstates are trivial and analytically given as \a Si S% S3 S M). The additional 
quantum number a refers to a set of intermediate spin quantum numbers which appears when 
coupling single spins of the same sublattice to the corresponding sublattice spins Si and further 
on coupling the sublattice spins to the total spin S. With regard to the set of intermediate spin 
quantum numbers a the number of all possible ways of constructing states characterized by fixed 
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values of the sublattice and total spin quantum numbers determine the degeneracy of energy levels 
in the spectrum of the rotational-band Hamiltonian. 

In the case of s = 1 the lowest band in Fig. [2]is given by states \a Si S2 S3 S M) with sublattice 
spin quantum numbers S\ = S2 = S3 — S max = 4-1 = 4 while the second band is given by a 
deviation of one sublattice spin of 1, i.e. \a (S'max — 1) <Smax S max S M) and permutations thereof. 
The other bands can then be constructed by introducing additional deviations of the sublattice 
spin quantum numbers. The energy spectrum with s = 3/2 shown in Fig. [2] can be constructed 
accordingly. 

Additionally several energy levels in Fig. [2] are colored. This coloring refers to so-called super- 
bands. A super-band consists of those rotational bands for which the sum of sublattice spin quantum 
numbers is the same. In contrast to bipartite systems (see Ref. [[6]) the spectrum of the rotational- 
band Hamiltonian for the cuboctahedron is much denser at low energies and only the first three 
super-bands are well separated from the others. 

The classical ground state of the system plays the key role within the approximate diagonal- 
ization. The better the quantum mechanical system can be approximated by a classical picture 
of the ground state the more effectively the approximate diagonalization works, i.e. the faster the 
approximate energy eigenvalues converge towards the exact values. From a purely classical point 
of view the ground state of a cuboctahedron exhibits a three-sublattice structure and is infinitely 
degenerate since there coexist coplanar and non-coplanar vector orientations [ [35] • Here it should 
be emphasized that within the approximate diagonalization, i.e. within a quantum mechanical 
treatment, coplanarity of the classical ground state does not play a role, but the coloring of the 
classical ground state does, since there is a direct impact on the set of intermediate spin quantum 
numbers a in the rotational-band states | a S\ S2 S3 S M) which are taken into account as basis 
states in the approximate diagonalization. 



r(q=0)-family M-family 




Figure 3. Families of the classical ground state of the cuboctahedron. Coloring of the Y{q = 0)- 
family (upper left). Coloring of the three equivalent states of the M-family (upper right and 
bottom). The drawing is a schematic planar projection. The labels correspond to y: yellow, r: 
red, and b: blue. 

The classical ground state of the cuboctahedron exhibits 24 colorings of the spins which can be 
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- by group theoretical considerations decomposed into two families that are invariant under oper- 
ation of the full point-group symmetry Oh [ HH] ■ Following Ref . [ [115] these families will be denoted 
as T(q = 0)-family and M-family. It has also been shown that those irreducible representations 
of Oh which form the T(q = 0)- and M-families are found in the low-lying part of the spectrum 
of the quantum cuboctahedron with half-integer spins s = 1/2,3/2,5/2 [[19]. Figure [3] shows the 
colorings of the different classical ground state families of the cuboctahedron. In a classical picture 
each color refers to a sublattice with all spins pointing into the same direction. The angle between 
classical spins belonging to different sublattices is 120°. 

In order to calculate the approximate spectrum of the system one is now left with the con- 
struction of basis states of the form \aS1S2S3SM), i.e. quasi-classical states. Therefore spins 
belonging to the same sublattice have to be coupled to yield the total sublattice spins Si, S2 and 
S3. Afterwards these total sublattice spins are coupled to the total spin S. The underlying cou- 
pling scheme is given by a classical ground state, i.e. a coloring from Fig. [3] and incorporated in 
the quantum number a. In the following the resulting basis states will be labelled with respect to 
the classical reference state. To this end one has to distinguish between basis states of the form 
1 7 Si S 2 S 3 SM) r , \ni Si S 2 S3 SM}m i: |a*2 Si S2 S3 SM)m 2 and |/i 3 Si S 2 S 3 SM)m 3 - The nota- 
tion of the set of intermediate quantum numbers and the subscript of the states now directly point 
to the classical ground state colorings, i.e. the underlying coupling scheme. It is important to note, 
that each of these four basis sets spans the same Hilbert space Tt{S, M). 

In the case of the M- family three different colorings exist. When diagonalizing the Hamiltonian 
approximately while additionally using point-group symmetries one has to restrict to those sym- 
metry groups where the symmetry operations do not alter the sublattice structure. A symmetry 
operation has no impact on the sublattice structure if the corresponding spin permutation results 
in recoloring of spins where all spins of a given sublattice maintain the same color, i.e. subscript (r: 
red, y: yellow, b: blue). For example an operation which leads to a cyclic permutation of colorings 
like 

r — > y — > b — > r (8) 

has no impact on the sublattice structure. 

While the sublattice structure of the T(q = 0)-family is left invariant under all symmetry 
operations of Oh, the sublattice structure (i.e. coloring) of the M -family is not. In order to restore 
the sublattice invariance of the classical ground state belonging to the M-family the basis states 
for the approximate diagonalization can be constructed from a usually overcomplete set M of basis 
states by an orthonormalization procedure. This set consists of rotational-band eigenstates from 
each coloring of the M-family, i.e. is given by 

M = { |/ii S 1 S 2 S 3 S M)\ h } + { |/i 2 Si S 2 S 3 S M)\ l2 } + { |^3 Si S 2 S 3 S M)\ h } , (9) 

while the index i is taking the values 1, . . . , n re d- Here n ro( j reflects the overall number of states 
contained in the incorporated rotational bands which is independent of the choice of the coupling 
scheme. 

Since the underlying coupling scheme is different regarding the chosen coloring these states 
have to be converted into an - in general - arbitrary reference scheme A before calculations can be 
performed. A transition between states belonging to different coupling schemes can be calculated 
using general recoupling coefficients [[33l[34] . For the transition of a basis state S M)m ± referring 
to a coloring Mi of the M-family into the scheme A one yields 

\aSM) A = A(aSM\ f M 1 SM) Ml \aSM) A , (10) 

OC,Sl,S2,S3 

where the summation indices indicate that the summation is running over all valid combinations 
of values for the intermediate spin quantum numbers a as well as for the sublattice spin quantum 
numbers. The transitions between states of colorings M 2 and M3 into the scheme A can be obtained 
analogously. 
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Figure 4. Classification by irreducible representations of Oh of the quasi-classical states with 
lowest energies of the T(q — 0)-family and M-family within S — and S = 1 subspaces. The 
given energies refer to a diagonalization of the Heisenberg Hamiltonian exclusively within these 
quasi-classical states in the system with s — 3/2. 



Figure 0] shows the classification of the quasi-classical states with lowest energies within 5 = 
and S = 1 subspaces according to irreducible representations of Oh- It corresponds to an ap- 
proximate diagonalization of ffneisenberg using only a reduced basis set of states from the lowest 
rotational band in the system with s = 3/2. Obviously the classification of the quasi-classical states 
is independent of the single spin quantum number s. Solely the energies arc changed when approx- 
imately diagonalizing the Hamiltonian of a system with different single spin quantum numbers. 
Based on a classical ground state belonging to the T(q — 0)-family the states belong to the irre- 
ducible representations A\ g (S = 0) and A 2g , E g (S=l). Taking the M-family as a starting point 
for the construction of basis states the states belong to Ai g , E g (S=0) and Ai g , E g , T 2u (S=l) 
where the T2 M -state has apart from its intrinsic degeneration, i.e. the dimension of the irreducible 
representation T 2ui a- twofold multiplicity. 

The difference in the geometric properties of the lowest states of the chosen ground-state family, 
expressed by their decomposition into different irreducible representations, should directly lead 
to a convergence behaviour which depends on the choice of the underlying coloring. Since the 
approximate basis set contains energetically ordered states (from lowest to higher), low- lying states 
are expected to converge faster with growing number of basis states used for the approximate 
diagonalization [[6]. 

3.1. Convergence of individual colorings 

In Fig.[5]thc approximate eigenvalues of a cuboctahedron with s = 1 and s = 3/2 in dependence 
on the number of incorporated bands within the S = subspaces are shown. As underlying classical 
ground states for each system the T(q = 0)-family as well as one member of the M-family have 
been chosen. In every subfigure the last column refers to results from a complete diagonalization 
within Ti.(S = 0). Since the T(q = 0)-family is invariant under Oh, the full point-group symmetry of 
a cuboctahedron (Oh) was used in order to classify the states. As mentioned before the individual 
members of the M-family are not invariant under Oh but under D 2 . Thus, in this case the D 2 
point-group symmetry was applied. 

Looking at Fig. [5] it becomes apparent that the convergence is directly dependent on the choice 
of the underlying classical ground states. In agreement with theoretical expectations states which 
correspond to low-lying eigenstates of the rotational-band model converge faster than higher-lying 
states. Additionally, by comparing the spectra of s = 1 and s = 3/2 it can be seen that the ground 
state energy converges more rapidly with increasing single spin quantum number s. Nevertheless, 
apart from some states which exhibit a quite regular and fast convergence behaviour the conver- 
gence is - overall - rather poor. Especially the so-called low-lying singlelQ, which appears in the 

x Here: the lowest E 9 -state in the system with s = 3/2 based on a classical ground state coloring of the F(q = 0)- 
family. 
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Figure 5. Low-lying energy spectra of a cuboctahedron with s = 1 (upper row) and s = 3/2 
(lower row) within S — subspaces using an increasing number of bands (states). The underlying 
classical ground states belong to the V(q = 0)-family (left) and to one coloring of the M-family 
(right) while Oh (left) and D2 (right) point-group symmetries were used in order to classify the 
states. The dashed boxes in the lower left figure indicate regions, i.e super-bands, in which the 
energy of the ground state is considerably lowered. The numbers refer to the number of magnons 
of the super-bands. 



case of half-integer spins s below the first triplet contributions, converges very slowly. This means 
that it is poorly approximated by one or a few eigenstates of the rotational-band Hamiltonian, but 
instead is a superposition of very many basis states. 

One significant property which becomes evident when tracing single energy levels is that the 
convergence is stepwise - at least in the low-lying part of the spectra. Regions in which the ground 
state is considerably lowered are excmplarily marked with dashed boxes in the case of a cubocta- 
hedron s = 3/2 assuming a classical reference state of the T(q = 0)-family. The numbers below the 
boxes refer to the number of magnons existing in the states of those rotational bands which are 
incorporated within the marked regions. Obviously the ground state energy is lowered whenever 
super-bands with an even number of magnons are incorporated in the approximate diagonalization. 
The observation of a stepwise convergence leads to a helpful additional approximation given by an 
approximate selection rule which is discussed below (see Sec. l3.4|l . It should be mentioned here that 
qualitatively the same convergence behaviour can be observed within subspaces of S > 0. How- 
ever, additional information cannot be extracted from graphical visualizations of the convergence 
in these subspaces, thus they will not be presented here. 

3.2. Convergence of combined colorings 

In order to improve the convergence of the energy levels in comparison to the behaviour shown 
in Fig.[S]the decomposition of the quasi-classical states in Fig. Indirectly leads to the starting point. 
Since the quasi-classical states, that result from different colorings, can be classified according to 
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irreducible representations of 2Ai g and E g , it is a straightforward task to use linear combinations 
of rotational-band eigenstates of the T(q = 0)-family as well as of the M- family as basis states for 
the approximate diagonalization. The set of basis states results from an extension of M in Eq. ([9]) 
to 

G = {\fn Si S 2 S 3 SM) l Ml } + {\n 2 Si S 2 S 3 SM)l h } 
+ { |M3 Si S 2 S 3 S M)\ h } + { | 7 Si S 2 S 3 S M)^} 

and subsequent orthonormalization of the incorporated states. 
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Figure 6. Energy difference to the ground state of selected low-lying energy levels of the cuboc- 
tahedron s = 1 within S = subspace in dependence on the number of incorporated rotational 
bands. The reduced set of basis states is taken from rotational-band states of the T(q = 0)-family 
only (left) and from linear combinations of states of all four colorings (right). The classification 
is according to Oh point-group symmetry. The dotted red line refers to the exact difference to 
the i? 9 -state. 
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In Fig. [S] the energy difference between selected low-lying energy levels and the ground state 
of the system is displayed for a cuboctahedron s = 1 within 5 = subspace. The difference is 
shown in dependence on the number of incorporated rotational bands. The last column refers in 
both subfigures to the exact values, taken from a complete diagonalization. The dotted red line 
indicates the energy difference to the £^-state. 

In the first case (left subfigure) only rotational-band states of the T(q = 0)-family are taken 
into account for the approximate diagonalization. In the second case (right subfigure) linear com- 
binations of states of all four colorings are used. The first noticeable difference is that when taking 
basis states of all colorings into account the low-lying levels start and remain in close proximity 
to their final (true) value. This is especially obvious when comparing the convergence of the low- 
est Aig-state in both subfigures. The second difference is that when using linear combinations of 
states of all four colorings the convergence is smoother and more rapidly. This can be traced back 
to the fact that the different colorings contribute low-lying states of different irreducible repre- 
sentations which in the other families would only be available as high-lying states. Therefore, it 
is advantageous to use fewer bands but basis states of all classical colorings for the approximate 
diagonalization. 

A short explanation might be appropriate in order to understand why the energy differences 
sometimes increase when taking more states into account. This is due to the fact that the ground 
state in such cases converges more rapidly than the excited states. Although all approximate 
eigenvalues improve on an absolute scale the differences get worse for a moment. 



3.3. Results 

In this section we like to present how good thermodynamic observables can be approximated. 
Figure [7] shows the specific heat C(T,B) (left) and the magnetic susceptibility dM/dB (right) for 
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zero field B — of a cuboctahedron s = 3/2. The energy spectrum of this system can be completely 
calculated using irreducible tensor operator technique in combination with a D2 point-group sym- 
metry [[35]. In Fig. [7] the exactly calculated specific heat and the susceptibility are compared with 
results from an approximate diagonalization using only states of the lowest rotational band in each 
coloring in order to set up the basis. We also show how these observables look like if evaluated 
with only the lowest rotational band (L-band) of Hrb (red colored energy levels in Fig. [21). 
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Figure 7. Comparison of approximate and exact specific heat C(T, B) (left) and zero-field suscep- 
tibility dM/dB (right) for a cuboctahedron a = 3/2. The L-band refers to magnetic properties 
evaluated with the lowest band of the rotational-band Hamiltonian. For the approximate di- 
agonalization only the states of the lowest rotational band but of each coloring are taken into 
account. 



As one can see the low-temperature specific heat is very sensitive to the structure of low- 
lying levels. In the rotational-band model the lowest band (L-band) is gaped and describes a 
rotor. Therefore, the specific heat is suppressed at small temperatures and at higher temperatures 
approaches 3/2 ke. The approximate diagonalization, although only taking the lowest (L-band) 
states of all four coloring into account, achieves already an improvement for low-temperatures. 
The little Schottky-peak is the result of low- lying singlets and rearranging triplets. At higher 
temperatures this approximation displays the same behavior as the pure L-band which is expected 
since only L-band like states are incorporated. 

In contrast to the specific heat which directly reflects the density of states within a certain 
energy interval the magnetic susceptibility reflects the density of magnetic states. To that end it 
is not surprising that the magnetic susceptibilities calculated from the approximate spectra do 
not differ considerably Since the exact spectrum also possesses a rotational-band like behaviour 
at the lower edge the contributions from the approximate diagonalization and from the L-band 
reproduce the steep rise of the susceptibility for low temperatures very well. Obviously, the exact 
thcrmodynamical properties cannot be reproduced properly in the case of increasing temperature 
because the approximate spectra only contain a fraction of the energy levels of the full Hilbert 
space H. 



3.4. Approximate selection rule 

As mentioned before (Sec. 13. ip the stepwise convergence behaviour leads to an approximate 
selection rule. Using this selection rule the computational effort when setting up the Hamilton 
matrices can be further reduced. In Fig. [5] it was shown that when diagonalizing in the T(q = 0)- 
family, regions can be marked in which single energy levels are affected by taking into account 
additional rotational-band states. The same can be done for the M-family. As it was already 
mentioned, the marked regions coincide with the incorporation of rotational-band states belonging 
to the same super-band. 

Looking at Fig. [3] the approximate energy eigenvalues of the ground state are lowered whenever 
states I7 S\ S*2 5*3 S M) are additionally included which belong to a super-band with an even number 
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of magnons. Since the convergence behaviour is state-sensitive the energy of the lowest states 
belonging to the 3-dimensional irreducible representation T\ u , for example, is considerably lowered 
whenever super-bands with an odd number of magnons are included. This relation between the 
incorporation of super-bands and the immediate affection on the energies of certain approximate 
eigenstates can also be observed if the classical ground state belongs to the M -family. 

As a result of the aforementioned observations a simple rule can be conjectured for low-lying 
energy levels |5i, 6 5*2,6 5*3,6 S M). The matrix elements (5i j0 S 2 , a 5 3 , a S M\H\Si,b S 2 ,b 5 3 ,& S M) con- 
necting these states with other states \S\^ a S2, a S3. a S M) are one magnitude (or more) bigger than 
other matrix elements if 

\n a -n b \ mod 2 = 0, (12) 

where n a = Xa=i 5i, ma x — Si=i &i,a an d n b (similarly evaluated), denote the number of magnons 
of the super-bands the states are belonging to. The Hamilton matrices approximately split up 
according to Eq. (Tl2^) . Thus, simultaneously computation time is saved and the dimensionality of 
the problem is reduced. 

4. Summary 

In this paper an approximate diagonalization scheme was proposed and used in order to deter- 
mine the energy spectra of geometrically frustrated spin systems. As an example the cuboctahedron 
was discussed. It was shown that the convergence is clearly dependent on the coloring of the un- 
derlying classical ground state, i.e on the coupling scheme. Therefore, the approximation can be 
improved by using an adapted set of basis states that originates from rotational-band states of all 
possible sublattice colorings. Furthermore, it was shown that the approximate diagonalization is 
useful in order to study the low-temperature thermodynamics of geometrically frustrated systems, 
even in its simplest form when only states from one rotational band within each coloring are taken 
into account. 

Although the necessary calculations are rather involved, especially for huge frustrated systems, 
the approximate diagonalization can be a valuable method. Numerically the evaluation of recou- 
pling coefficients constitutes the strongest challenge, i.e. calculations are rather limited by runtime 
than by available memory. However, recent developments show that highly parallelized program 
code or public resource computing can help to overcome this barrier. 
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